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Abatraet— We  present  results  of  experiments  and  numerical  simulations  designed  to  reveal  the  presence  of  * 
spontaneous  oscillations  arts.ng  from  negative  differential  mobility  effects  in  gallium  arsenide  field  effect  Iran-  v 
sistors  The  measurements  include  d  c  and  pulsed  current/voltage  vs  temperature  charactenzation.  sampling  scope  u* 
measurements,  spectral  analysis  to  40GHi  and  observation  of  light  emission  The  simulation  is  a  time  dependent 
large  signal  transient  analysis  arising  from  a  fully  two-dimensional  solution  of  the  self-consistent  potential  and 
charge  within  the  device 
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The  gallium  arsenide  field  effect  transistor  has  been 
around  a  long  tune  (1.2)  showing  promise  of  delivering 
high  power  while  sustaining  low  noise  levels.  GaAs  is. 
however,  subject  to  hot  electron  effects  that  manifest 
themselves  in  the  presence  of  highly  nonuniform  charge 
layers  or  "domains"  And  depending  on  device  size  and 
bias  values,  these  domains  are  either  trapped  near  the 
gate  or  dram  contacts(3).  or  propagate  between  them. 

The  presence  of  trapped  high  field  domains  is  expected 
to  profoundly  affect  FET  performance  and  is  currently 
subject  to  much  study  (4-6).  However,  with  the  excep¬ 
tion  of  early  empirical  studies  [2. 7)  and  recent  numerical 
studies(3,8. 10).  the  dynamic  properties  of  the  propagat¬ 
ing  domain,  the  conditions  for  their  existence,  and  their 
influence  on  device  performance  has  until  recently III- 
14)  been  largely  ignored  The  renewed  interest  arises 
from  a  need  to  determine  the  extent  to  which  the  occur¬ 
rence  of  spontaneous  oscillations  in  GaAs  FETs[ll-l4] 
are  a  consequence  of  negative  differential  mobility 
(NDM)  In  the  following  we  discuss  the  influence  of 
trapped  and  propagating  domains  on  FET  operation.  We 
also  present  evidence  for  the  highly  probable  conclusion 
that  the  propagating  domain  is  the  origin  of  our  obser¬ 
vation  of  spontaneous  oscillations  in  GaAs  FETs. 

The  study  reported  here  is  in  two  parts:  (i)  numerical 
simulation  and  (ii)  experiment  The  numerical  studies 
reveal  the  presence  of  two  qualitatively  different  trapped 
domains  and  one  traveling  domain  The  first  trapped 
domain  forms  within  the  conducting  channel  and  at  the 
drain  edge  of  the  gate  contact.  [4J6]  Its  formation  is  a 
direct  consequence  of  velocity  limitation  in  GaAs  and  its 
presence  is  responsible  for  current  saturation  at  values 
of  drain  voltage  significantly  below  that  predicted  by 
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ShodUey[l5).  The  second  trapped  domain  is  qualita¬ 
tively  different  and  forms  in  FETs  capable  of 
sustaining  a  current  instability  arising  from  (NDM).  This 
domain  forms  near  the  drain  contact,  at  high  drain  and 
low  gale  bias  levels  and  is  coincident  with  the  cessation 
of  oscillatory  activity  in  GaAs  FETs.  (The  extent  to 
which  it  is  removed  from  the  gate  region  depends  on  the 
gate  to  drain  spacing  and  the  doping  concentration  of  the 
epitaxial  region.)  The  domain  is  qualitatively  similar  to 
the  high  anode  field  configuration  in  two  terminal 
devices,  where  it  is  also  responsible,  in  many  cases,  for 
the  cessation  of  oscillation  (16).  Each  of  the  low  and  high 
drain  bias  trapped  domains  may  pulsate  in  time  and  so 
enhance  the  noise  existent  in  the  FET.  In  the  case  of  the 
traveling  domain,  it.  as  for  some  transferred  electron 
logic  devices  (TELD)  is  launched  at  the  drain  aide  of  the 
gate  contact,  propagates  to  the  drain  contact  and  then 
recycles  successively  While  the  traveling  domain  need 
not  be  "stable"  in  the  sense  ’of  two  terminal  NDM 
devices  (16).  its  threshold  and  quenching  conditions  are 
similar  to  that  for  two  terminal  devices  (3, 8). 

The  experimental  studies  reported  below  consist  of 
pulsed  and  d.c.  measurements  performed  as  a  function  of 
temperature  and  under  different  circuit  conditions.  We 
have  determined  the  spectral  content  of  the  oscillations 
when  they  occurred,  and  for  pulsed  measurements  have 
sampled  the  current  and  voltage  at  different  points  within 
the  pulse  The  devices  studied  were  empirically 
separated  into  two  groups  as  determined  by  the  ratio 

K- 

_ Gate  voltage  at  cutoff  ( ■  V„) _ 

Drain  voltage  at  the  onset  of  current  saturation  ( ■  y*)‘ 

(1) 

Devices  with  K  greater  than  unity  sustained  current 
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oscillations.  those  with  k  approumatels  equal  to  or  less 
than  units  were  electrically  quiet  l  or  those  cases  with  k 
creator  than  unitv  oscillation'  began  bevond  the  knee  of 
current  saturation  In  each  case  the  oscillations  ceased  at 
high  \ allies  of  drain  and  (tale  has  In  most  of  the 
electric alh  unstable  devices  the  onset  of  the  oscillation 
is  marked  h\  a  drop  in  current  In  all  cases  the  instability 
appoirs  lt>  be  initiates I  when  th <  aieraite  heU  unJcr  the 
v'sitc  Cs'ntaet  first  rcsiches  the  NPM  lhresh<>U  hehl  rails' 
This  result  holds  true  under  different  ambient  tem¬ 
perature  conditions  An  estimate  of  the  carrier  velocitv  in 
the  channel  between  the  cate  and  drain  contacts  in 
dicates  that  prior  to  the  oscillation  the  electrons  arc 
ususills  sin  ft  me  sit  spcs'sls  in  excess  of  the  s,iturati<’n  slnft 
iclocits  When  the  oscillations  ceased  with  increasing! 
cate  bias  the  carrier  vekvitv  dropped  to  values 
significantlv  below  the  saturation  drift  vekvitv  bach  of 
the  above  results  are  consistent  with  an  instabilitv  arising 
from  the  presence  of  negative  differential  mobilitv 

:  nr  vii  m  i  vssirK  vrios 

The  classification  parameter  k  i  see  cgn  1 1  emphasizes 
differences  between  the  l«.i  and  classical  I  FT.  where 
for  the  latter 1 1 V.  and  I are  equal  For  Ga As.  the 
verv  earliest  experiment'll.  2. '|  showed  k  v  .using  from 
unitv  to  values  significantlv  higher  than  unitv  i('urrent 
dev  ices |l'|  have  ni't  altered  this  situation*  In  those 
e\perimenis| 1. 21  when  k  •  I.  no  instabilhv  was  repor 
ted  When  k  was  at  the  high  end  of  the  scale,  the  studies 
in  Ref'  |2|  and  PI  showed  that  the  device  no  longer 
exhibited  the  classical  characteristics,  but  instead  sit' 
t.uned  spontaneous  oscillations  for  a  range  of  bias 
value' 

The  result'  me  computer  simulation  corroborate  the 
empirical  cla-  situation  Rut.  in  addition  relate  stabilitv  to 
the  material  parameter  r„  the  electron  saturated  drift 
vekvitv  SpcciAcallv .  the  simulations  show  that  a  nee*'- 
sarv  condition  for  a  current  instabilitv  is  that  the  current  at 
saturation  either  exceed  or  be  within  a  narrow  band  about 
the  value 

i.  =  OsF.i  <:» 

In  eqn  (2>.  G,.  is  the  open  channel  conductance 
i=  Area  i  L  '  Resistixitxi).  E.  -  r.  a.  L  is  the  source  to 
drain  separation,  and  u  the  low  held  mobilitv  of  the 
semiconductor  In  two  terminal  device  language.  F  is  the 
sustaining  held|lb|  Summarizing  If  IDDS  denotes  the 
zero  gate  bias  current  at  saturation  then 

IDSS  s- 1.  Stable  current-voltage  characteristic' 

IDSS >  /. :  Instable  current-voltage  characteristics.  i?l 

The  abv've  conclusion  is  essentiallv  similar  to  that  of 
Yamaguchi  et  al  (8] 

The  result'  discussed  above  summarize  the  general 
features  of  the  Ga  As  FFT  The  space  charge  distribution 
within  the  FET  is  responsible  for  these  properties  and 
will  be  discussed  in  the  next  section 
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Mil)  Introiluction 

The  following  analvsts  is  for  devices  in  the  configura¬ 
tion  of  big  Mat  where  the  source  and  drain  contacts  are  al 
parallel  ends  ol  the  device  This  configuration  eliminates 
geometrical  effects  due  to  coplanat  placement  of  the 
source  and  dram  contacts  I'oplanar  contact  calculations 
have  been  performed  for  the  configuration  shown  in  lig 
Itbi.  which  also  include  the  presence  of  a  substrate  And 
the  results  of  these  studies  'how  infection  into  the  sub 
strate.  with  ihe  instabilities  still  remaining  Ml  of  the 
c  ilculations  reported  beiow  are  for  elements  with  doping 
levels  nominallx  equal  to  IO,\'ir.\  a  value  approx.  I  * 
to  2  orders  of  magnitude  below  that  of  presentli  fahri 
cated  devices  I  his  choice  was  made  because  Ihe  cost  of 
computing  ihe  contributions  of  transiting  domains  in 
creases  monotonicallv  with  increased  doping  level | IH| 
However,  a  number  of  representative  calculations  have 
been  performed  at  higher  doping  levels  with  similar 
results  We  note  that  the  results  of  Yamaguchi  el  al  |8| 
were  at  *  x  10"'  cm1  and  where  there  is  overlap  with  our 
studies  there  is  agreement  Also,  the  mam  results  of  the 
experiment  confirm  the  conclusion  developed  at  these 
lower  doping  levels 

kb)  Current  roltaec  characlen sties  anil  space  charge 
Jistnbulion  for  bb'T  with  k  I 

The  space  charge  distributions  for  a  k  I  device  are 
displaced  in  bigs  24  The  inset  of  each  figure  shows  its 
current -voltage  characteristic  with  the  lower  case  letters 
keved  to  the  space  charge  distribution  The  current  is 
normalized  to  the  v  .due 

lr  =  GiiF.,.1..  (4) 

where  Fr  i'  the  electric  field  it  peak  vekvitv  And  the 
voltage  is  normalized  to 

Vr  -  F,  l  .  <5t 

The  bold  curve  is  the  electron  drift  vekvitv  electric  field 
relation  scaled  to  current  and  voltage  The  closed  circles 
denote  computed  point'  and  each  point  is  a  stable  time 
independent  point 

For  the  sp.ue  charge  distribution  the  r-axis  represent' 
the  longitudinal  spatial  dependence  The  vertical  charge 
distribution  i'  along  the  i  -axis  The  source  contact  is  m 
the  v-r  plane  at  x  =  0:  the  dram  contact  is  at .«  =  /  The 
gate  contact  is  in  the  i plane  at  the  indicated  points. 

At  low  values  of  drain  bias  and  zero  gate  bias  (Fig  2a). 
there  is  a  region  of  charge  depletion  under  the  gate 
contact  The  charge  region  spreads  somewhat  hevond  the 
gate  contact  boundaries  (enclosed  rectangular  parallele¬ 
piped).  but  prettv  much  the  space  charge  distribution  is 
what  one  would  expect  from  classical  arguments|l5]  An 
increase  in  drain  voltage  to  the  knee  of  the  currerd 
voltage  relation  (Fig  2bi  leads  to  a  domain,  representing 
the  first  significant  departure  from  the  classical  analvsts 
The  domain  consists  of  an  accumulation  laser,  which 
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t  ig.  I.  Device-circuit  configuration  used  in  the  FET  simulations,  (a)  I  arallel  source  and  drain  contacts,  (h)  Coplanar 
source  and  drain  contacts.  Substrate  is  included 


forms  within  the  boundaries  of  the  gale  region  and  is 
followed  downstream  by  a  depletion  layer.  Domain  for¬ 
mation  is  a  consequence  of  current  continuity  and  velo¬ 
city' limitation  and  occurs  when  the  current  density  under 
the  gate  exceeds 

h  =  AU'p  (6) 

whereiy,  is  the  maximum  carrier  velocity.  Further  in¬ 
creases  in  drain  bias  result  in  an  increase  in  the  amount 

of  charge  accumulated  and  a  broadening  of  the  domain 
(Fig.  2c).  We  point  out  that  while  the  presence  of  a 
domain  is  necessary  for  a  current  instability,  it  is  not 
sufficient.  By  analogy  with  two  terminal  devices  the 
current  density  downstream  from  the  dipole  must  exceed 
the  sustaining  current,  7„  where 

7,  =  N0ev,  (7) 

(see  also  eqn  2).  This  criteria  is  approximately  satisfied 
for  the  FET  but  is  not  as  stringent,  It  suffices  to  state 
that  the  current  density  for  the  Ac  ==  f  devices  is  generally 
too  low  to  sustain  an  instability. 

The  calculations  of  Fig.  3  are  for  finite  values  of  gate 
bias.  At  low  drain  bias  levels  the  depletion  layer  extends 
well  beyond  the  gate  boundaries.  We  see  a  region  of 
local  charge  accumulation,  although  clearly  there  is  a  net 
depletion  of  carriers.  However,  increasing  the  drain  bias 
results  in  a  strong  region  of  charge  accumulation,  fol¬ 


lowed  by  significant  depletion  until  the  drain  contact  is 
reached  where  the  boundary  conditions  require  a 
significant  increase  in  the  charge  density. 

Two  extreme  conditions  are  displayed  in  Fig.  4.  Here 
at  low  values  of  drain  bias  and  a  value  of  gate  bias  near 
cutoff  we  see  the  presence  of  a  depletion  layer  broadly 
surrounding  the  gate  region.  At  zero  gate  bias  levels  but 
very  high  drain  bias  levels  we  see  the  presence  of  a 
region  of  charge  accumulation  extending  from  the  gate 
region  to  the  drain  contact. 

We  summarize  the  above  results  by  stating  that  for  a 
collection  of  ostensibly  classical  current-voltage  curves 
the  space  charge  distribution  is  far  different  in  many 
circumstances  from  that  developed  originally  for  the 
FET(  15).  Thus  while  the  current-voltage  relation  does 
not  provide  us  with  detailed  information  about  the  dis¬ 
tribution  of  space  charge  for  the  K  **  I  device,  we  expect 
that  the  form  of  the  distribution  will  affect  the  equivalent 
circuit  element  representations  of  the  device  [4,5], 

3(c).  Current  voltage  characteristics  and  space  charge 
distributions  for  FETs  with  K  >  1 

For  a  given  doping  level  and  source  drain  spacing,  the 
principle  difference  between  devices  with  K  >  1  and 
K  =  I  is  that  the  former  have  wider  channels,  therefore, 
lower  open  channel  resistances,  and  for  a  given  bias 
draw  more  current.  For  these  devices  the  current  density 
is  usually  high  enough  to  sustain  an  instability.  In  our 
simulations  the  onset  of  the  instability  marked  the  onset 
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pi_  2  SDace  charKe  configuration  as  a  function  of  bias.  Profiles  correspond  to  indicated  points  on  the  current  voltage 

characteristics  Solid  line  represents  GaAs  velocity  vs  electric  field  relation  scaled  to  current  and  voltage  Devtce 
characteristics,  ooi  confiBVlrat.on  is  also  shown,  with  dimensions  in  microns. 


of  current  saturation.  In  our  experiments,  the  instability 
occurs  after  the  device  enters  saturation,  The  simula¬ 
tions  for  the  wider  channel  device  are  displayed  in  Figs. 

5  and  6.  . .  .  . 

Figure  5(a)  shows  the  charge  distribution  at  a  bias  level 
just  below  that  necessary  to  launch  a  traveling  domain. 
Within  the  gate  boundaries  there  is  charge  depletion.  But 
towards  the  bottom  of  the  channel  and  just  before  the 
boundary  there  is  a  trace  of  an  accumulation  layer- 
implying  that  the  current  density  is  highest  there.  Coup¬ 
ling  this  result  to  the  cun-ent  value  in  the  inset  to  Fig.  5, 
we  have  the  result  that  within  the  gate  boundary  the 


current  density  is  approximately  equal  to  JP,  while  within 
the  gate-to-drain  region  the  current  density  exceeds  J,. 

The  above  conditions  are  sufficient  to  launch  a  domain 
in  two-terminal  devices,  and  indeed  here  a  slight  increase 
in  drain  bias  results  in  the  launching  of  a  domain  (Fig. 
5b).  In  this  case  further  increases  in  bias  result  in  the 
charge  accumulation  layer  spreading  toward  the  bottom 
of  the  channel  (within  the  gate  boundaries)  and  also 
toward  the  drain  contact.  The  domain  detaches  from  the 
gate  region  and  periodic  transit  time  oscillations  occur. 
Figure  5(b)  shows  a  "snapshot"  of  the  domain  at  one 
instant  of  time  during  propagation.  (Recycling  domains 
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Fig  3  As  in  Fig  2.  hut  for  different  bias  values 


are  illustrated  in  Fig  5.  Ref  (3).)  The  propagating  and  goes  to  zero  and  the  system  relaxes  to  the  value  denoted 

recycling  domain  can  persist  for  a  rather  wide  range  of  by  the  closed  circle  When  the  domain  is  propagating  the 

drain  bias  values  The  amplitude  of  these  resulting  transit  current  oscillates  along  the  load  line.  At  a  bias  of  1.6  V'p 

time  oscillations  are  bias  dependent  as  illustrated  in  Fig.  the  current  oscillation  amplitude  is  somewhere  around 

7.  20 9c  of  its  d.c.  value.  Further  increases  in  drain  bias 

Figure  7  contains  plots  of  drain  and  source  current  vs  result  in  an  increase  in  the  amplitude  of  the  oscillation, 

drain  potential  at  different  times,  with  time  eliminated  The  average  current  associated  with  these  oscillations  is 

between  these  In  Fig  7(a)  we  plot  drain  current  vs  drain  denoted  by  the  x's.  We  note  that  the  average-current  - 

voltage  for  four  different  values  of  drain  bias  The  bold  voltage  relation  for  the  range  in  which  there  are  tramit 

lines  denote  calculated  results  the  dashed  lines  are  an  time  oscillations,  exhibits  a  region  of  negative  differential 

estimate  of  the  current  vs  voltage:  the  skew  lines  are  the  conductivity  (NDC).  The  NE  it,  for  these  calculations 

drain  d.c.  load  lines  In  each  calculation  displacement  dynamic  in  origin.  The  second  part  of  Fig.  7  shows  the 

current  contributions  result  in  transient  drain  current  source  current  vs  drain  voltage  for  the  same  calculation 

values  in  excess  of  their  ste.idv  state  values.  This  is  as  that  of  7(a).  Here,  in  Fig.  7(b)  we  see  the  presence  of 

apparent  for  the  bias  levels  of  1.4  and  I  6  V'p.  where  after  superlinear  current-voltage  relations  at  low  drain  bias 

steady  state  has  been  reached  the  displacement  current  levels,  followed  by  looping  when  domains  are  present. 
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Both  of  these  effects  occur  because  of  large  transients 
through  the  gate  contact. 

With  respect  to  transit  time  oscillations  we  point  out 
that  no  experiment  has  vet  been  designed  that  pros  ides 
unequivocal  evidence  for  their  existence  in  commonly 
fabricated  GaAs  RET"  In  most  common  designs  the  gate 
to  dram  separation  is  somewhere  between  one-half  and 
five  microns,  and  for  these  values  it  is  not  clear  that 
conditions  exist  for  the  propagation  of  "stable" 
domains|l6)  The  result  is  that  the  properties  of  these 
propagating  domains  are  readily  influenced  by  the 
external  circuit  and  can  in  some  cases  be  controlled  bv  it 
An  example  of  this  complex  dipole-circuit  interaction 
exists  in  two  terminal  devices  (e  g  the  quenched  multiple 
dipole  model  Ifijl  There  is  also  the  possibility  that  in 
place  of  the  negative  differential  mobility  induced  oscil¬ 


lations  in  GaAs  FETs  there  are  feedback  circuit  oscil¬ 
lations  that  are  independent  of  any  hot  electron  effects. 
However,  in  the  experiments  reported  below  we  show 
that  the  evidence  points  to  the  oscillations  as  being 
negative  differential  mobility  induced  Two  additional 
space  charge  profiles  are  important  in  this  regard  and  are 
discussed  below 

The  two  profiles  of  interest  are  coincident  with  the 
suppression  of  transit  lime  oscillations.  In  one  case  a 
significant  increase  in  drain  bias  alone  suppresses  the 
instability.  In  the  second  case  an  increase  in  the  mag¬ 
nitude  of  the  gate  bias  is  sufficient  to  suppress  it  The 
space  charge  profile  associated  with  the  high  dram  bias 
suppression  of  the  instability  is  displayed  in  Rig  fi(a)  We 
see  that  an  accumulation  of  charge  forms  downstream 
from  the  gate  boundaries  and  extends  towards  the  drain 
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contact.  While  this  layer  generally  is  not  a  static  one  hut 
instead  pulsates  in  time  there  are  some  broad  features 
associated  with  it  For  example,  the  profile  forms  at 
rallies  of  current  density  in  excess  of  J,  and  at  hiph 
rallies  of  drain  potential.  Most  of  the  potential  drop  is 
across  the  gate-to-drain  region  ( Ref. (3|.  Fig.  7),  where  we 
ma\  expect  the  electric  fields  to  he  well  within  the 
saturated  drift  velocilx  region  Thus  the  carriers  within 
the  gate-to-drain  region  are  traveling  at  their  saturated 
drift  velocity  values.  This  fact  and  current  continuity 
with  (he  current  density  exceeding  ./,  are  responsible  for 
the  region  of  charge  accumulation.  By  way  of  analogy 
this  solution  is  qualitatively  similar  to  two  terminal 
device  solutions  which  show  the  presence  of  anode 
adjacent  domains)  16).  Further,  in  three  terminal  devices 
with  long  gate-to-drain  spacings  and  high  doping  con¬ 


centrations  the  anode  adjacent  domain  has  been 
experimentally  probed|ld|. 

The  space  charge  profile  associated  with  the  high  gate 
bias  suppression  of  the  instability  is  displayed  in  Fig. 
htb)  This  shows  a  snapshot  of  the  space  charge  layer 
while  undergoing  a  damped  oscillation  This  layer  of 
charge  also  pulsates  in  time,  but  in  contrast  to  the  space 
charge  laser  of  Fig  Wa).  this  one  forms  at  current 
density  rallies  that  are  less  than  ),.  Indeed  for  this  case 
most  of  the  potential  drop  falls  under  the  gate  contact 
and  extends  some  distance  downstream  from  the  edge  of 
the  gate  boundaries  |1|  (Fig.  8)  As  long  as  the  space 
charge  layer  does  not  reach  the  drain  contact  the  electric 
field  downstream  from  the  edge  of  the  nonuniform  space 
charge  region  will  be  less  than  the  "two-terminal  sus¬ 
taining  field"  |  lf»|  FV  =  i  .Ip  This  situation  is  analogous 
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to  the  two  terminal  high  cathode  boundary  field 
profiles  1 16]. 

3(d).  Summary 

To  summarize,  there  are  several  key  features  asso¬ 
ciated  with  negative  differential  mobility  induced  in¬ 
stabilities  in  field  effect  transistors.  These  are: 

(1)  The  instability  occurs  in  devices  with  K  greater 
than  unity. 

(2)  The  instability  occurs  at  or  beyond  the  “knee"  of 
the  cunent-voltage  relation 

(3)  A  necessary  condition  for  an  instability  is  that  the 
current  density  downstream  from  the  gate  boundary  be 
greater  than  or  approximately  equal  to  J„  or 
equivalently,  that  the  drift  velocity  cither  exceeds  v, .  or 
is  approximate^  equal  tc  r. 

(4)  Saturation  in  current  occurs  when  (he  average  field 
within  the  conducting  channel  under  the  g-ite  contact  is 


approximately  equal  to  the  threshold  field  for  negative 
differential  mobility. 

(5)  The  instability  can  be  suppressed  by  imposing  a 
gate  bias  sufficiently  high  to  reduce  the  current  density 
downstream  from  the  gate  boundan  lo  a  value 
significantly  less  than  J„ 

(6)  The  instability  can  be  suppressed  by  imposing  a 
sufficiently  high  value  of  drain  bias. 

Let  us  now  turn  our  attention  to  how  these  instabilities 
manifest  themselves  in  actual  devices. 

5  EXPERIMENTS  WITH  GAIJJl'M  ARSENIDE 
FIELD  EFFECT  TRANSISTORS 

5(a).  Experiments  with  K  >  I 

Figures  II  and  9  display,  respectively,  temperature 
dependent,  d.c.  and  pulsed  drain  current  vs  drain  voltage 
characteristics  for  a  GaAs  FET.  The  source,  drain  and 
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!2  kV/cnt,  and  we  have  reduced  these  results  by  approx. 
10%  for  higher  fields.  The  closed  circles  in  Fig.  10(a)  are 
for  the  d.c.  measurements;  the  x's  represent  pulsed 
measurements  ( I  p  sec  pulse  lengths,  1%  duty  cycle).  The 
pulsed  results  generally  surround  the  solid  line  and  are 
higher  than  the  d.c.  results.  In  general  over  the  temperature 
range  exceeding  150  K  the  pulsed  measurements  yield  a 
velocity  exceeding  i\. 

(4)  The  average  field  under  the  gate  boundary  is  esti¬ 
mated  from  the  relation 


Ek  =  V,JLk  (9) 

where  VlK  is  the  longi'udinal  potential  drop  under  the 
gate  contact  and  is  obtained  from  an  equation  that 
ignores  contact  resistance  contributions: 


OOO  0.25  0.50  0.75  1.00  1.25  1.50  1.75 

vONAtN/'vP 

Fig  7.  Instantaneous  current  vs  voltage  for  the  device  of  Figs.  5 
and  6.  (al  (Train  current  vs  drain  voltage  (b)  Source  current  vs 
drain  voltage. 

gate  contacts  are  coplanar  with  the  source  and  drain 
contacts  separated  by  approx  8.5  ^m.  The  gate  is  more 
or  less  centrally  placed  with  a  length,  L„.  equal  to 
3.0  /cm.  The  gate  width,  Z.;  is  250  jcm.  The  epitaxial  layer 
thickness,  H,  was  3000  ±  500°  A  with  a  nominal  doping 
level  of  I0,1/cm\  The  mobility  measured  from  adjacent 
samples  varied  from  3000  to  4000cm2/V-s.  The  cutoff 
voltage  for  this  device  was  8  9  V.  This  device  sustained 
current  oscillations  as  detected  by  a  spectrum  analyzer. 
All  six  criteria  listed  in  the  last  section  for  NDM  induced 
instabilities  were  satisfied.  We  discuss  each  separately 
below.  Oscillation  characteristics 

(1)  K=  V'„,/V'</<a<8.9  +  0.8V)/’.5  =  3.9,  where  the 
cutoff  voltage  includes  as  assumed  built  in  potential  of 
approx.  0.8  V. 

(2)  At  each  value  of  ambient  temperature,  the  in¬ 
stability  occurred  beyond  the  knee  of  the  current  voltage 
characteristic. 

(3)  The  velocity  was  computed  from  the  relation 


r  -  l'lNueHL.).  (8) 


V,„  =  Vj.  -IRJL-IJIL  (10) 

where  is  the  drain  voltage  and  /  the  drain  current. 
Ro{=  I/O,,)  is  the  open  channel  resistance  and  may  be 
obtained  from  the  measured  drain  current  voltage 
characteristics  of  the  device  via  the  relation 

Ro  =  RA  V„)( I  - 1  +  V„)l  V,„]l,;) /(l  “  f  I  - ^] ( V* 

+  VH)IV,„)m)  (11) 

where  we  have  adopted  the  viewpoint  of  eqns  (2)  and  (3) 
of  Grebene[22).  V,„  in  eqn  (II)  includes  the  built  in 
potential.  The  above  expression  for  R0  is  expected  to  be 
valid  at  low  values  of  drain  bias.  Assuming  a  built  in 
potential  of  0.8  V  and  reading  Rd  at  V,  =  0  from  the  d.c. 
characteristics  we  compute  R0.  The  values  are  listed  in 
Table  1,  where  we  have  ignored  the  temperature  varia¬ 
tions  in  V*  and  V„„  We  see  that  the  average  field  under 
the  gate  is  relatively  insensitive  to  temperature,  as  is  the 
threshold  field  for  negative  differential  mobility [20],  An 
average  field  of  4.2kV/cm  is  consistent  with  a  mobility 
for  IOl7/cm'  material  of  4200cm2/V-s. 

(5)  As  shown  in  Fig.  10(b)  the  carrier  drift  velocity  at 
the  cessation  of  oscillations  is  below  v,. 

(f>)  All  oscillations  cease  at  high  drain  bias  levels. 

The  next  set  of  data  are  displayed  in  Fig.  1 1  where  we 
show  d.c.  and  sampling  scope  data  (taken  25  nsec  into  a 
50  nsec  pulse,  with  low  duty  cycle).  The  GaAs  FET  had  a 
source  to  drain  spacing  of  3.5  pm  and  a  gate  length  of 
I  pm.  The  cutoff  voltage  for  this  device  was  3.0V  and 
the  mobility  was  quoted  at  around  4500cm2/V-s.  With 
regard  to  the  first  two  oscillatory  criteria,  K  ~  4.8  and  the 
instability  occurs  just  beyond  the  knee  of  the  current  vs 
voltage  characteristic.  The  peak  velocity  just  prior  to  the 
instability  at  57  ma  is  computed  from  the  equation 


And  the  values  prior  to  an  instability  are  displayed  as  a 
function  of  ambient  temperature  in  Fig.  10(a)  where  the 
error  bars  represent  uncertainties  in  the  value  of  H.  In 
Fig  10(a)  the  solid  line  corresponds  to  an  estimate  of  the 
saturated  drift  velocity  from  Ruch  and  Fawcett[20|  and 
Freeman  and  Hobson[2l|.  This  data  is  available  only  to 


v  =  plRoL  (12) 

with  Rn  =  9.7 11  (from  pulsed  data)  and  p  =  4500cm2/V- 
s,  The  result  is  e -0.68  x  I07  cm/sec.  which  falls  ap¬ 
proximately  within  the  range  of  values  of  Fig.  10  An 
average  threshold  field  of  4  kV/cm  for  a  one  micron  gate 
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ached  at  a  drain  voltage  of  approx.  0.8  V;  and  it  is 
that  we  see  substantial  preinstabihty  current 
dening.  The  oscillations  which  are  initiated  ut  about 
I  are  not  suppressed  until  a  sufficiently  low  value  of 
ent  density  is  attained.  The  characteristic  of  the 
ibility  is  suggestive  of  hot  electron  contributions 
er  than  that  discussed  in  the  last  section  (see  also 
(13])  However,  consistent  with  the  sixth  criteria  the 
nations  are  suppressed  at  high  values  of  drain  bias. 
,  is  displayed  in  Fig.  12  where  we  show  sampling 
r*  measurements  at'  increasing  dram  bias  for  the 
ice  subjected  a  l.^sec  pulse.  We  see  the  initiation 


of  oscillations  at  low  drain  bias  levels  with  subsequent 
suppression  at  high  drain  bias  levels. 


cl.  Experiment  with  K  ~  I 

The  final  set  of  measurements  is  for  a  device  wit 
fiproximately  equal  to  unity.  These  are  pu  se 
leasurements.  The  device  has  a  cutoff  voltage  o  p- 
rox  I  8  V  and  no  oscillations  were  detected  on  a  spe  - 
■urn  analyzer  over  a  range  of  40  GHz  The  device  had  a 
uoted  doping  concentration  of  1  *  10  /cm  and  an  epi¬ 
cal  thickness  of  0.5  gm.  These  values  were  no. 
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Fig  <t  Temperature  dependent  pulsed  data  for  the  deuce  of  Fig.  X  For  these  measurements  a  <  f !  resistor  sepa  at.  d  the 

source  from  ground 


consistent  with  pinchoff  values  using  the  equation 

=  «'N'oW'7(2<).  (13) 

With  a  quoted  value  of  / ...  =  360 gem.  eqn  (8)  yields  a  room 
temperature  range  of  values  of  t  .  associated  with  the 
product  ,V„H. 

r  =  4.7  x  |0IH/M,tf.  (14) 

From  eqn  (14)  with  /V„W  varying  from  IO'J/cm2  to  1.5  x 
10';/cm‘,  f  varies  from  about  0.5  x  107  cm/sec  to  0.3  > 
I07  cm/sec.  All  vttlues  are  too  low  for  an  instability  to 
occur. 

4(c).  Discussion 

The  experiments  discussed  above  provide  evidence 
that  the  GaAs  FFT  is  capable  of  sustaining  negative 
differential  mobility  induced  current  instabilities.  These 
instabilities  occur  when  the  average  field  under  the  gate 
contact  exceeds  the  NDM  threshold  field  and  when  the 
current  density  at  the  onset  of  saturation  is  ap¬ 
proximately  equal  to  or  exceeds  the  saturation  drift 
current  density  for  GaAs.  There  are  several  additional 
experimental  aspects  that  should  be  highlighted.  The  d.c, 
measurements  of  Fig.  8  were  characterized  by  a  drop  in 
current  at  the  instability  threshold:  the  pulsed  measure¬ 
ments  also  show  a  drop  in  current.  The  d.c.  measurements 
of  Fig.  1 1  do  not  show  a  dropback  in  current;  the  sampling 
scope  measurements  do.  There  is  hysteresis  in  the  d.c. 
measurements  of  Fig.  8  with  the  instability  persisting  for 
increasing  and  decreasing  values  of  drain  bias.  However 
on  the  decrease  the  instability  persisted  at  drain  bias 
levels  somewhat  below  the  instability  threshold.  This  is  a 


characteristic  of  two  terminal  transit  linv  uc.ia-.i If*!. 
With  regard  to  the  sampling  scope  measurements  of  Fig. 
II.  we  see  some  spreading  of  the  current  \s  voltage  prior 
to  the  instability.  This  may  correspond  in  p;rt  to  a 
pulsating  trapped  domain  near  the  gate  contact,  as  dis¬ 
cussed  in  the  earlier  sections. 

The  oscillation  properties  reported  here  v.o  c  detected 
using  sampling  oscilloscopes  and  spectrum  analy/eis. 
Spectral  analysis  was  used  on  the  devices  ol  Figs  8.  9 
and  13.  Sampling  scope  measurements  were  for  the 
device  shown  in  Figs.  II  and  12  and  other  devices  not 
reported  here.  At  first  the  spectrum  analyzer  was  set  up 
to  detect  frequencies  to  18GHz.  In  this  set  up.  frequen 
cies  of  1.2,  1.6,  2.5  and  4.2 GHz  were  detected.  These 
oscillations  were  bias  dependent  hut  ceased  at  drain  bias 
levels  in  excess  of  6  V.  The  oscillations  ceased  it  sub¬ 
stantially  higher  gate  bias  levels.  Under  the  assumption 
that  there  were  additional  frequencies  in  excess  of 
18  GHz  that  were  going  undetecled.  and  that  by  analogy 
to  two  terminal  devices  where  the  introduction  of  a 
inductor  reduces  the  circuit  frequency,  we  introduced  a 
JnH  inductor  in  series  with  the  source  loop.  The  low- 
frequency  oscillations  disappeared  and  in  their  place 
were  instabilities  at  13,  13.6  and  !4.2GHz.  More  detailed 
spectral  analysis  measurements  were  subsequently  made 
with  the  spectrum  analyzer  open  to  frequencies  to 
40 GHz.  The  pulsed  measurements  were  then  analyzed, 
with  the  following  details.  Oscillations  began  at  3.5  V  and 
persisted  to  ju.'  below  4.25  V.  Oscillations  were  obser¬ 
ved  at  6.5,  13,  23,  30,  36  and  40GHz.  Higher  frequencies 
could  not  be  detected.  The  oscillations  were  strongest  at 
4.1V  where  the  6,5  and  4AGHz  signals  dominated  At 
4.25  V  on  the  drain  the  oscillations  ceased,  at  least  up  to 
40GHz.  In  the  same  measuring  circuit  and  with  the  gate 
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Fig  10  Solid  curve  is  an  estimate  of  high  field  saturated  velocity  vs  amhent  temperature  tat  Velocity,  from  eqn  8.  vs 
ambient  tempera'-  -'■.prioi  to  oscillations,  (hi  Velocity .  from  eqn  8,  vs  ambient  temperature,  at  cessation  of  oscillations 

with  application  of  a  high  eate  bias. 


wide  open  visible  light  was  detected  at  6.4  V  op  the 
drain.  This  was  not  accompanied  by  any  oscillations 
whose  frequency  was  below  40 GHz.  An  initial  reduction 
in  potential  on  the  gate  contact  resulted  in  similar  behavior 
but  at  -0.8  V  on  the  gate, all  oscillatory  activity  ceased, 
including  the  appearance  of  light.  The  electrical  behavior 
of  this  device  was  measured  at  a  variety  of  temperatures, 
with  qualitatively  similar  behavior. 

The  oscillation  frequencies  reported  above  bear  no 
clear  relation  to  the  transit-time  oscillation  frequency.  As 
indicated  in  Section  3.  the  gate  drain  separation  may  be 
too  short  to  guarantee  stable  domain  propagation  and  so 
the  oscillatory  properties,  including  the  frequency  are 
likely  to  be,  determined  by  a  complex  device-circuit 
interaction.  There  is  also  the  possibility  that  the  oscil¬ 
lations  we  are  observing  are  a  consequence  of  feedback 


between  the  gate  and  drain  loop.  While  this  situation 
cannot  be  ignored  the  evidence  in  terms  of  the  six 
instability  criteria  given  in  the  last  section,  including  the 
suppression  of  oscillations  at  high  drain  bias  levels,  make 
this  possibility  small 

In  the  above  paragraphs  we  reported  the  observation 
of  light.  Light  was  also  observed  unuer  d  c.  conditions. 
The  details  are  as  follows.  Under  d  c.  conditions  light 
was  seen  at  the  beginning  of  the  current  instability, 
became  dimmer  as  the  drain  bias  increased  but  remained 
while  the  device  oscillated.  The  light  was  stronger  at 
drain  bias  levels  sufficient  to  eliminate  the  oscillations. 
At  moderate  gate  bias  levels  faint  light  remained,  but  a! 
sufficiently  high  gate  bias  levels  large  enough  to  eliminate 
the  oscillations,  ail  light  activity  ceased.  The  light  for  all 
piactical  purposes  was  white  The  observations  of  light  is 
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Fig  1 1 . (a)  d.c. and;(b)  pulsed  sampling  scope  measurements  of  the  current  voltage  characteristics  for  a  short  gate  FET 
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Fig  13.  Pulsed  current  vs  voltage  data  for  a  device  that  did  not  exhibit  any  oscillations 


compelling  evidence  for  the  pttscnce  of  regions  of  high 
electric  field  and  excess  earners.  At  low  drain  bias  levels, 
the  dimness  of  the  light  may  be  accounted  for  by  the  fact 
that  the  time  averaged  electric  field  at  a  given  point 
within  the  sample  is  less  for  propagating  domains  then 
for  the  trapped  domains.  The  latter  occurs  at  high  drain 
bias  levels.  The  absence  of  light  at  high  gate  bias  levels 
may  possibly  be  due  to  the  broad  region  of  charge 
depletion  in  the  vicinity  of  the  gate  boundary  The 
arguments  for  white  light  are  somewhat  uncertain.  We 
recall,  however,  that  the  separation  between  the  central 
portion  of  the  conduction  band  and  the  valence  band  is 
within  the  infrared,  while  the  separation  between  the 
satellite  conduction  and  valence  bands  is  within  the 
ultraviolet  The  entire  visible  spectra  is  within  these  two 
bounds 

We  note  that  Mimura  et  al.  [23]  on  IO'’/cm'  GaAs 
FETs  also  observed  v  isible  white  light  at  high  drain  bias 
levels  For  devices  with  highly  doped  n '  drain  contacts 
light  was  observed  at  the  n'n  interface;  for  devices 
without  an  n'  drain  contact  radiation  was  observed  at 
the  drain  metalization  edge  Mimura  et  a/  (23).  attribute 
the  radiation  to  impact  ionization  within  a  domain  trap 
ped  at  the  drain  contact  The  intensity  of  this  light  was 
ooserved  to  increase  with  drain  bias  and  decrease  as  the 
gate  voltage  was  made  negative  These  results  are  con¬ 
sistent  with  our  findings 

5  SOW.  WOOCIJJV.  (  ONNIMJUTtOSS 

The  results  of  our  analysis  are  consistent  with  the 
identity 

(15) 


where  is  the  longitudinal  electric  field  in  the  source 
to  gate  region  In  devices  with  K  greater  than  unity, 
will  generally  be  at  least  equal  to  E, ;  and  so  for  this  case 
a  lower  bound  for  IDSS  exists: 

IDSS  *  Rn  S  E,L.  (16) 

For  example,  with  the  device  reported  in  Figs.  8  and  9, 
the  room  temperature  values  of  the  l.h.s.  of  eqn  (16) 
(from  Table  I)  is  equal  to  1.9V  The  rh.s.  side  with 
£,  =  2.2  kV/cm  is  equal  to  1.87  V.  For  the  device  whose 
properties  are  displayed  in  Figs  II  and  12,  IDSS  x  Ra  = 
0.52  V.  while  E,L  =  0.77  V.  By  way  of  comparison 
Shockley's  study  (24)  teaches  that 

IDSSxR 0-VJ3  (17) 

while  more  recently  Shur (25)  have  concluded  that 

IDSS  x  R„  =  E.L,.  (18) 

Equation  (18)  yields  results  significantly  lower  than 
observations,  and  less  than  eqn  ( 16). 

t  CONdl  SIONS 

In  conclusion,  through  analysis  and  experiment  we 
have  presented  evidence  for  the  presence  of  negative 
differential  mobility  induced  current  instabilities  in  GaAs 
FETs.  We  have  listed  six  criteria  that  identify  these 
NDM  instabilities,  paramount  among  these  is  the  fact 
that  the  longitudinal  electric  field  under  the  gate  must  be 
greater  than  the  threshold  field  for  negative  differential 
mobility  while  the  current  density  must  exceed  the  sus- 
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(aiding  current.  The  fact  that  the  instabilities  in  some 
cases  are  delayed  until  higher  values  of  bias  are  reached 
implies  that  the  field  must  exceed  the  NDM  threshold 
field  for  the  instability.  We  have  also  presented  cor¬ 
roborating  experimental  and  theoretical  evidence,  con¬ 
sistent  with  the  conclusions  of  Yamaguchi  et  al.  [8],  that 
devices  whose  current  densities,  within  the  source  to  gate 
region,  are  less  than  J,  will  not  go  unstable.  Those  whose 
current  densities  exceed  J,  will  sustain  instabilities. 
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